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Can long-term periodic variability and jet helicity in 
3C 120 be explained by jet precession? 
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ABSTRACT 

Optical variability of 3C 120 is discussed in the framework of jet precession. Specifi- 
cally, we assume that the observed long-term periodic variability is produced by the 
emission from an underlying jet with a time-dependent boosting factor driven by 
precession. The differences in the apparent velocities of the different superluminal 
components in the milliarcsecond jet can also be explained by the precession model as 
being related to changes in the viewing angle. The evolution of the jet components has 
been used to determine the parameters of the precession model, which also reproduce 
the helical structure seen at large scales. Among the possible mechanisms that could 
produce jet precession, we consider that 3C 120 harbours a super-massive black hole 
binary system in its nuclear region and that torques induced by misalignment between 
the accretion disc and the orbital plane of the secondary black hole are responsible for 
this precession; we estimated upper and lower limits for the black holes masses and 
their mean separation. 

Key words: galaxies: active - galaxies: individual: (3C 120) - galaxies: jets - radio 
continuum: galaxies 



1 INTRODUCTION 



3C120 (z=0.033; iBaldwin et"all Il980t) . also known as II 
Zw 14 and PKS 0430-1-052, is usually classified as a 
Seyfert 1 galaxy, although its morphology in the optical 
band is not as simple as that of a typical galaxy of this 
class. Indeed, photometric and spectroscopic studies seem 
to indicate that 3C 120 either pas sed, or it is still pass- 
in g through a mer ger process (e.g., ISoubevran et alJll989l: 



I^ortl^^^ 



Hiorth et al 



Il995l) The residual I-band image obtained by 
after subtraction of the stellar con- 
tribution, showed a complex structure formed by several 
condens ations, probably associa ted to active star forming 
regions iSoubevran et all Il989^ , and an elongated struc- 
ture which coincides with the kilo-parsec rad io jet detected 
at 5 GHz JWalker. Walker fc Bensonl Il988^ . In fact, this 
large scale jet is the extension of a sub-parsec scale jet, 
which remains relativistic up to d istances of about 100 kpc. 
JWalker. Benson fc Unwinlll987a^ . 

Several superluminal radio components have been de- 
tected in the jet (e.g ., Gomez ct al. 1998, 2000; Walk er et alj 
I2OOII : IComez et alj|200lf) . with different velocities and po- 
sition angles, b esides a stationary core smaller th an 54 pas 
(0.025/1"^ pc) iComez. Marscher fc Alberdilll999f) . There is 
also evidence of the existence of trailing shocks in the jet 
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jGomez et al.ll200lD . that is, features that do not originate 
in the jet inlet. These features could be related to pinch- 
mode jet-body instabilities produced by the propagation of 
the superluminal com ponents, as shown r ecently by numer - 
ical simulations (e.g.. lAgudo et alJlioOll : lAlov et al.ll2003f) . 
The observed correlation between dips in the X-ray emis- 
sion and ejection of superluminal components has been in- 
terpreted as a consequence of the connection between jet 
origin a nd the accretion disc, such as in the case of micro- 
quasars iMarscher et alj|2002f) . 

30 120 prese nts variability in all bands and in d i fferent 
time-scales (e.g., |EBste in_etal ]|l972l: lHalDernlll985t IWebbI 



ll990l:IShukla fc Stoneill996l : IZdziarski fc Grandi2001l) . The 

complex variability found in the histo rical _B-band light 
curve was decomposed by IWebbI il990fl in three different 
components: a linear time decrease of its magnitude due to 
the diminution of the accretion rate, a long-term variabil- 
ity with a period of 12.4-yr produced by thermal or viscous 
instabilities in the accretion disc, and short-term variations 
associated to magnetic eruptions in the magnetized disc. 

In this paper, the reported 12.4-yr variability is in- 
terpreted as periodic boosting of the radiation emitted by 
the underlying jet, caused by jet precession. This model 
also explains the differences in superluminal velocities 
and position angles of the different components, assuming 
that they represent the direction of the jet inlet at the 
epoch in which the components were formed. Further- 
more, the complex jet structure of 3C 120 at large scales 
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is also studied in the framework of jet precession. Jet 
precession has been claimed by several authors in order to 
explain the radio st ructure o f several quasars, BL Lacs and 
radio -R alaxies (e.g..lGower fc Hutchin^ s 1982: Gower et alJ 
19821: iGower fc Hutchingj Il984l: Iroos fc Meuerj 



AbrahamI 12001 ISthling et ali 120031: ICaproni fc AbrahamI 
200^ )~"suggesting that jet precession is not so uncommon 
phenomenon in the Universe. We shall adopt h = 0.7 and 
go = 0.5 throughout the paper. 



2 JET PRECESSION 

2.1 Observational evidences 

HiRh-resolutio n obse rv ations between 5 and 4 3 GHz 
JWalker et aP Il982l: IWalker. Benson fc Un^ Il987b l: 
Gomez. Marscher fc Alberdil Il999l: iGomez et all I200C 



Fomalont et al.l l2000l: iHoman etal] l200ll: IWalker et alJ 
20^F ~g'ive the core-component distance r and the posi- 
tion angle on the plane of sky rj for each feature in the 
sub-parsec scale jet. Using data at different epochs, the 
apparent proper motion n is calculated and the apparent 
velocity /3app in units of light speed c is obtained from: 



/5app — 



goz+(go-l)[(l + 2go2)^/"-l] 



(1) 



where h = _ffo/100, Hq is the Hubble constant in units of 
km s~^ Mpc~^, go is the deceleration parameter and z is 
the redshift. The ejection epoch to of each component is 
obtained by back-extrapolation of their linear motions. The 
kinematic parameters of the different superluminal features 
are presented in Table 1. 

We have labelled jet components as 'K' followed by a 
number related to the epoch in which they were formed 
('1' for the oldest one). We also present the labels given 
in earlier works in the second column of Table 1. Except for 
KIO, we kept strictly the previous identifications. Consider- 
ing the uncertainties, the listed parame ters are compatible 
with those found in the literature (e.g.. iGomez et aLlll998l 
I2OOII: I Walker et alJl200lh . 

We can note in Table 1 that the different jet compo- 
nents were ejected with difi'erent velocities and position an- 
gles. A possible interpretation for this behaviour is in terms 
of a precessing jet model: precession changes the orientation 
of the jet inlet in relation to the line of sight, so that the 
direction in which the components are ejected, as well as 
their apparent velocities become a function of time. For the 
present discussion, it is not relevant whether the jet com- 
ponents are assumed to be plasmons or shocks since we are 
only interested in the kinematic aspects. 

At lower resolution, VLBI observations at 1.7 GHz have 
reve aled that the jet str ucture of 3C 120 is extremely com- 
plex JWalker et alJ200J) : it presents sub-structures in scales 
of tenths of a parsec that move superluminally, a possi- 
ble stationary component located at an angular distance 
of about 81 mas from the core and a jet aperture that is 
larger in the southern direction, specially after about 180 
mas, where there is also a decrease in the app arent velocity. 

Those characteristics were interpreted bv lWalker et alJ 



i2001f) as an indicative of the presence of a helical pattern 
in the jet. However, they can also be understood in terms of 
the precession model that explains the sub-parsec behavior 
of the superluminal jet. 

2.2 Precession model 

We will derive the instantaneous appearance of the jet, as- 
suming that it is the result of the combination of plasma 
elements ejected in different epochs, with different angles in 
relation to the line of sight. Let us consider a plasma ele- 
ment ejected at time to with velocity C/3 (c is the light speed) 
in the comoving reference frame. This element will have a 
velocity c/3app in the observer's reference frame given by: 



/3app — 



±P8m[4>{to)] 
l±/3cos[0(to)] 



(2) 



where (p is the angle between the moving direction of the 
element and the line of sight. The signs and '-' refer 
to the jet and the counterjet respectively. As the counter- 
jot has not been detected for scales smaller than 100 kpc 
(Walker. Benson & Unwin 1987a), we will consider only the 
jet hereafter. 

Due to jet precession, and rj are functions of time t 
given by: 



(j}{t) — arcsin 

ri{t) — arctan 
with 



^x{t)^+y{ty 



x(t) 



^{t) — A{t) cos 770 — B{t) sin 770 



y{t) = A{t) sin 770 4- B{t) cos 770 



and 



A{t) — cos n sin 00 + sin i} cos <po sin uit 



Bit) — sinficostjt 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



where uj is the precession angular velocity, is the semi- 
aperture angle of the precession cone, 0o is the angle between 
the precession cone axis and the line of sight and 770 the 
projected angle of the cone axis on the plane of the sky. 

In the determination of the precession parameters, we 
used all jet components listed in Table 1 as model con- 
straints. As data were obtained at different frequencies, 
shifts in the core-component distances due to opacity ef- 
fects may become quantitatively important in the calcula- 
tion of the proper motions. To perform opacity correction 
in the observational data, we used the formalism given in 
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Table 1. Kinematic parameters of the superluminal components of 3C 120. 



Component 


Literature'' 


to (yr) 


fi (mas 


■/yr) 




V (°) 


Kl 




1976.7 


± 


0.6 


3.01 


± 


0.33 


4.6 


± 


0.5 


-113 lb 2 


K2 




1977.6 


± 


0.6 


3.01 


± 


0.39 


4.6 


lb 


0.6 


-108 lb 4 


K3 


All] 


1978.8 


± 


0.5 


2.75 


lb 


0.33 


4.2 


± 


0.5 


-101 lb 9 


K4 


B[i] 


1980.3 


± 


0.5 


3.01 


it 


0.26 


4.6 


± 


0.4 


-99 ± 5 


K5 


c[i] 


1981.0 


± 


0.5 


2.48 


lb 


0.26 


3.8 


± 


0.4 


-95 ± 6 


K6 




1981.7 


± 


0.4 


2.35 


± 


0.20 


3.6 


± 


0.3 


-97 ± 6 


K7 


E[i] 


1982.6 


± 


0.4 


2.22 


± 


0.20 


3.4 


± 


0.3 


-101 lb 7 


K8 


F[il 


1983.3 


± 


0.5 


2.09 


± 


0.20 


3.2 


± 


0.3 


-102 lb 8 


K9 


G[il 


1984.3 


± 


0.5 


2.09 


± 


0.39 


3.2 


lb 


0.6 


-109 lb 8 


KIO 


H+I[il 


1985.4 


± 


0.8 


2.42 


± 


0.39 


3.7 


± 


0.6 


-119 lb 5 


Kll 




1986.1 


± 


0.6 


2.81 


± 


0.33 


4.3 


± 


0.5 


-120 ± 6 


K12 


JW 


1986.7 


± 


0.6 


2.35 


± 


0.33 


3.6 


± 


0.5 


-120 ± 7 


K13 


Kill 


1988.0 


± 


0.4 


2.81 


± 


0.39 


4.3 


± 


0.6 


-117 ± 7 


K14 


API 


1994.5 


± 


0.7 


2.22 


± 


0.26 


3.4 


± 


0.4 


-108 ± 3 


K15 


bPI, K1A/U1A131 


1994.9 


± 


0.7 


2.16 


± 


0.26 


3.3 


± 


0.4 


-106 ± 3 


K16 


CPl, KlB/UlB^l 


1995.2 


± 


0.6 


2.09 


± 


0.26 


3.2 


± 


0.4 


-111 ± 3 


K17 


£,[2,4]^ ^[2,4] 


1995.5 


± 


0.4 


2.16 


± 


0.20 


3.3 


± 


0.3 


-113 ± 3 


K18 


02121, gl21 


1996.3 


± 


0.5 


2.09 


± 


0.26 


3.2 


± 


0.4 


-119 ± 5 


K19 


H[2,4]^ 


1996.8 


± 


0.4 


2.09 


lb 


0.20 


3.2 


± 


0.3 


-120 ± 5 


K20 


j[2,41 j[2] 


1997.0 


± 


0.4 


2.16 


± 


0.26 


3.3 


± 


0.4 


-116 ± 3 


K21 


Kl2A][ k[5l 


1997.3 


± 


0.4 


2.22 


± 


0.20 


3.4 


± 


0.3 


-117 ± 5 


K22 


L12.41, 12[5l 


1997.5 


± 


0.4 


2.29 


± 


0.20 


3.5 


± 


0.3 


-124 ± 4 


K23 


ol+o2l51 


1998.2 


± 


0.4 


2.22 


± 


0.13 


3.4 


± 


0.2 


-122 ± 3 



Nomenclature in pre vious papers: [^yyalker et alj i200l|): PjQ^j^g^^ Marscher fc Alberd] OUS)); 
P lHoman et"ai] <200lD : ["^ iGomez etal\ i200(t : IS lGomez et al.l i200lD . 



iBlandford fc KonieJ l)l979|) . which uses the integrated syn- 
chrotron luminosity Lgyn, the ratio between upper and lower 
limits of the energy distribution in the relativistic jet par- 
ticles Fmax/rmin, the iutriusic jet aperture angle a con- 
stant parameter kc and the angle </!> between the jet direction 
and the line of sight. The relation between these quantities 
is present ed in Appendix A. 

As in lLobanovl il998h . we assumed kc — 1, Fmax/rmin = 
100, = 0.5° and Lsy„ = 8.4 x 10*^ erg s"^ However, in 
our model, the angle <j> is a. function of time and depends 
on the precession model parameters. Therefore, the correc- 
tions were determined iteratively together with the model. 
The final result was a small correction, with a mean value 
between 5 and 43 GHz of Arcoro ~ 0.147 mas, with upper 
and lower limits of 0.130 and 0.164 mas, respectively. 

To find the precession parameters we consider a period 
of 12.3 yr^, al most the sam e as the long-term variability 
period found bv lWebbI Jl990l) in the B-band light curve. In 
order to determine the best set of model parameters, we 
chose a 7 factor compatible with the velocity of the fastest 
component (~ 4.6^"^ c), given by -ymin = (1 + /^app)^^^ 
{'ymin ~ 6.7 for h — 0.7). After fixing a value for 7 close to 
its lower limit, we selected the parameters Q, (j>o and 770 that 
fitted the apparent velocities and position angles of the jet 
components. Then, we checked the behavior of /i/3app and 
r] as functions of time, assuming that the position angle of 
the jet component represents the jet position at the epoch 



^ In the source's reference frame, the precession period corre- 
sponds to 11.9 yr 



Table 2. Parameters of the precession model for the parsec-jet 
of 3C 120. 



P (yr)=' 7 


n n 


4>o (°) 


m n 


12.3 lb 0.3 6.8 lb 0.5 


1.5 lb 0.3 


4.8 ± 0.5 


-108 ± 4 



Measured in the framework fixed at the observer. 



when the component was formed. This procedure was done 
iteratively until a good fitting for the data was obtained. 

The precession parameters are given in Table 2, while 
the model fitting in the {rj, /i/3app), [t, 'i/^app) and {t, rj) 
planes is presented in Fig. 1. It is also possible to fit the 
data using larger values for 7 and other jet parameters (de- 
creasing fl and (po). As a consequence of that, not only the 
predicted position angles are much larger than those ob- 
served in the VLBI maps but also the time variations of 
the Doppler boosting factor become smaller than those nec- 
essary to explain the optical light curve (see Section 4 for 
further discussion). 



3 LARGE SCALE JET STRUCTURE OF 3C 120 

Taking into account that even though the jet components 
have ballistic trajectories on the plane of the sky, a snapshot 
in time of a precessing jet will show an helicoidal pattern. 
To compar e our results with th e large-scale jet structure de- 
scribed bv I Walker eEaDlllool, we calculated the apparent 
proper motion of a jet element from cf3a.pp[4>{t)] through 
equation (1). 
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Figure 2. Helicoidal pattern due to jet precession. Eacli one of 
the panels corresponds to snapshots of the precession helix in five 
distinct ep ochs, chosen as the same epochs of the 1.7 GHz ob- 
servations <Walker et'aI]|2001^ . The big circles superposed upon 
the helices represent the position of jet components which c an be 
related to the features L2-L9 found bv I Walker et al.1 120 



Figure 1. Precession model for 3C 120. The parameters are listed 
in Table 2. (a) - (c): Solutions of the model on the planes (rj, 
h/3app), (t, h/3app) and (t, rj). 

Using equations (2)-(8), we determined the right ascen- 
sion and dechnation offsets of the jet element in relation 
to the core (Aa and A5 respectively) in a given time tobs 
(tobs > to) through: 

Aa{tohs) = ^J■a{to) ■ (tohs - to) (9) 



A5(iob.) = Mi(io) • (iobs -to) (10) 

where fia and fj,s are respectively the apparent proper mo- 
tions in right ascension and declination, being related to /i 
and -q by: 

fia{to) = fJ,(to) ■ sin[r?(io)] (11) 



fJ.s{to) = m(*o) ■ cos[77(to)] (12) 

Assuming now a continuous jet and the precession pa- 
rameters listed in Table 2, we simulated the jet appear- 
ance at the five disti nct epochs for whic h 1.7 GHz obser- 
vations are available JWalker et al.ll200j) : t^bs = 1982.77, 
1984.26, 1989.85, 1994.44 and 1997.70, as shown in Fig. 2. 
Note that this approach is similar to that used to model the 
radio structur e of several quasars and radio-galaxies in pre- 
vious papers jGower fc Hutchingj|l982l: Ig ower et alJll982l : 
iGower fc Hutchingj|l984HRoos fc Meuerj|l987^ ■ 

Each point of the helices in Fig. 2 corresponds to a given 
plasma element, although in the case of a continuous jet, the 
helix driven by precession is continuous. The wavelength of 
the helices is basically related to the precession period, while 
the amplitude depends mainly on Q and (po. As this pattern 
is not stationary in time, since the precession helix is tied 
up with the jet movement, the configuration observed in a 
given time t will be seen again after an interval At, which 
corresponds to 12.3 yr in the case of 3C 120. 

The comparison of the predicted helicoidal patt e rns in 
Fig. 2 with the VLBI maps obtained bv lWalker et alJ (1200 J) 
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Table 3. Kinematic parameters of the superluminal components 
observed at 1.7 GHz in the large scale jet of 3C 120. 



Component 


*o (yr) 


M (mas/yr) 


/l/3app 


'7{°) 


L2 


1921.3 


2.24 


3.4 


-95 


L3 


1945.0 


2.49 


3.8 


-90 


L4 


1957.2 


2.52 


3.9 


-90 


L5 


1966.4 


2.94 


4.5 


-101 


L6 


1970.2 


2.32 


3.5 


-93 


L7 


1980.7 


2.76 


4.2 


-91 


L8 


1983.7 


2.10 


3.2 


-108 


L9 


1992.1 


2.88 


4.4 


-95 



show that the observed jet aperture is reasonably well de- 
scribed by jet precession in the inner parts of the jet. The 
maps reveal also the existence of discrete structures labelled 
as L2-L9, which could be superluminal. In order to compare 
these observations with our model, we assumed that L2-L9 
were ejected at different epochs and have ballistic motions. 
We present in Table 3 their kinematic parameters, which 
are used to calculated their right ascension and declination 
offsets for the five epochs given in Fig. 2 (their positions are 
marked by big circles) . Components L7 and L8 can be iden- 
tified respectively as the evolved components K5 and K8 in 
the parsec-scale jet, as presented in Table 1. 

A good agreement between the offsets shown in Fig. 2 
and the locations of L2-L9 in the 1.7 GHz maps of 3C 120 is 
found, suggesting that the simple approach assumed in this 
work (ballistic motion -\- precession) provide a reasonable de- 
scription for their kinematic behaviour in such scales. How- 
ever, this conclusion could be somehow misleading, since we 
can rule out neither the possibility that L2-L9 are formed 
by superposition of several unresolved components in the 
maps, nor that jet components have non-ballistic trajecto- 
ries in those scales. 

At larger distances from the core (| Aa |> 80 mas), we 
see that the full range of position angles provided by our 
precession model is not found in the observational data, in 
the sense that the amplitudes of the precession helices to- 
wards the southern direction are systematically larger than 
the observed jet aperture. A possible explanation for this 
behavior is the existence of an external medium which does 
not allow the jet propagation in some directions and could 
lead to the formation of a stationary component, as observed 
in the VLBI maps at Aa ~ —80 mas and A5 « — 17 mas. 

From these results, it is reasonable to verify in what 
conditions the quasi-stationary component could be formed. 
The interaction between the fluids at different velocities will 
produce shock waves propagating along the jet. In the case of 
strong shocks and assuming that the fluid can be described 
by a relativistic adiabatic equation of state (p oc N'^^^ , 
where p is the thermodynamic pressure), the jump condi- 
tion between the prc-shock and pos-shock regions is (e.g., 
iBlandford k McKca..l976. : .RomcrQ.199&'l : 



7^^ + (7ps7)^ 



(13) 



where Nj and 7 are respectively the proper particle density 
and the bulk Lorentz factor in the pre-shock region, while 
N^s and 7ps are the same quantities in the pos-shock region. 



If the observed feature is actually stationary, its proper 
motion must not produce displacements larger than the an- 
gular resolution of the maps, otherwise its motion would 
have been detected during this interval. From the beam-size 
in right ascension and declination (4 and 12.5 mas respec- 
tively) and the interval between the last and first observation 
(~ 14.93 yr), we estimated an upper limit for the apparent 
proper motion of about 0.27 mas/yr, which results in an 
apparent velocity smaller than 0.41c. 

To calculate the true velocity of the quasi-stationary 
knot, it is necessary to know its viewing angle. However, 
only its position angle (about —102°) is known and the pre- 
cession model allows two different solutions for the viewing 
angle: 3.3° and 6.2°. Using equation (2) with = 6.2°, we 
found that the upper limits for the true velocity and the as- 
sociated Lorentz factor for the quasi-stationary component 
are respectively 0.80c and 1.67, while for cj> — 3.3° the same 
quantities have upper limits of 0.88c and 2.11. 

With these limits and assuming that the quasi- 
stationary component is associated with the pos-shock re- 
gion, we found from equation (13) that the minimum particle 
density required to slow down the jet is 12A'j and lOA'j for 
the viewing angles 6.2° and 3.3°, respectively. For an exter- 
nal medium like the clouds of the Narrow Line Region, with 
typical number density of 10* cm^"^, and for a viewing angle 
of 6.2°, we obtain an upper limit of of 833 cm""' for the jet 
density, which is reasonable in terms of ordinary relativistic 
iets fe.g IWalker. Benson fc Unwinlll987al : lAltschuleJ[T98gl : 
iRomeroiligg^ 



4 THE OPTICAL BOOSTED EMISSION OF 
THE UNDERLYING JET 

The luminosity of 3C 120 observed in the optical band might 
be the result of the superposition of thermal contribution 
from the accretion disc with non-thermal contribution pro- 
duced, e.g., in the underlying jet and in the superluminal 
components. If we consider that the underlying jet is rep- 
resented by a relativistic continuous fluid, its fiux density 
measured in the observer's reference fra me S^iv) is related 
to that in the comoving frame S'j{u) by jLind fc BlandfordI 



2 + a 



(14) 



where a is the spectral index (5*1, oc u °') and 5 is the 
Doppler factor, defined as: 



5((^,7) = [7(1 cos </>)]- 



(15) 



Due to jet precession, 5, and consequently the observed 
emission from the underlying jet, become time dependent; 
the closer the jet is to the line of sight, the higher the boost- 
ing effect is. 

Depending on the intrinsic jet intensity, periodic 
boosting variations could produce detectable variability 
in the hght curve, as in the case of 3C 2 79 and OJ 287 
J Abraham fc Carraral IT99i: lAbrahamll2000t) . To rep r oduce 
this light curve and following the suggestion bv lWebbl(ll99(i) 
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Table 4. Model parameters used in equation (17) to describe 
long-term variability in the _B-band light curve of 3C 120. 



S[ (MJy) S'^ (nJy) 

MODEL A 0.19 -0.14 

MODEL B 0.16 0.0 



cases, the secular linea r decrease in the B-band brightness 
proposed originally by IWebbI il99Cf) seems to be unneces- 
sary, at least after 1965. We show in Fig. 3 the flux density 
calculated from model B, the B-band historical light curve 
of 3C 120 and the difference between them. 




1965 1975 1985 1995 

Epoch (years) 



2005 



Figure 3. Contribution of the underlying jet for the variability 
in the B-band light curve of 3C 120 in the case of Model B (see 
Table 4). The full circles are related to the observed flux density, 
while open squares are the residuals obtained after subtracting 
the contribution of the underlying jet (solid lines). The dashed 
line corresponds to residual equals to zero. 



of a possible secular linear variation, we modeled the con- 
tribution of the underlying jet assuming that its emission is 
described by equation (14), with S'^{v) given by: 



S[ [v, t) = S[ [v) + S'i (//) • {t - 1965) 



(16) 



The term 5(0,7)7(1 + z) in equation (16) is due to the 
transformation of the time measured in the comoving frame 
to the observer's reference frame. Combining equations (14) 
and (16), we have: 



■(1 + ^) 



2 + a 



(17) 



Observations in the B and / bands of the optical coun- 
terpart of the kilo -parsec radio jet led to as-i ~ 2.0 
l)Hiorth et alll995h : assuming that this value is also valid for 
the parsec-scale jet and using 5 from the precession model, 
we simulated the long-term periodic variability in two dif- 
ferent cases: Model A, with a secular time decrease of the 
intrinsic flux density, and Model B, in which the linear term 
was neglected (^2 ~ 0). The model parameters are given 
in Table 4, where S[{i') was chosen to keep Si(v,t) smaller 
than the minimum observed value at all epochs. For that 
reason the values of the parameters should be considered as 
upper limits. 

Boosting effects produce a substantial increase in the 
underlying jet flux density in both models. Since the residu- 
als obtained from the subtraction of the observed flux den- 
sity from the model predictions are very similar in the two 



5 A POSSIBLE SUPER-MASSIVE BLACK 

HOLE BINARY SYSTEM IN THE NUCLEAR 
REGION OF 3C 120 

Jet precession can be pr oduced by the Lense-Thirring effect 
iLense fc ThirriIl^ll918^ ■ in which precession is due to the 
misalignment between the angular momenta of the accre- 
tion disk and of a Kerr black hole. This effect was investi- 
gated bv ILuI il992h for several AGNs; for all of them, the 
period found was of the order of thousands of years. Using 
the same assumptions for the central object in 3C 120, whic h 
has a central mass of 3.4 X lO'^ Mq iPeterson et al.lll998l) . 
we obtain an even longer precession period. On the other 
hand, jet precession with periods of several years can be 
produced in super-massive black hole binary systems, when 
the secondary black hole has an orbit non-coplanar with 
the primary accretion disc, which induces torques in its in- 
ner parts (e.g., .Katz..l98Q , 1997 : .Romero et al.ll2000l) Thus, 
we will assume that the former scenario, in which jet inlet 
precession is induced in super-massive black hole binary sys- 
tem, is more suitable to 3C 120 and, from this assumption, 
the binary system parameters will be estimated. 

Let us consider that the primary and secondary black 
holes, with masses Mp and Ms respectively, are separated 
by a distance r^s- From Kepler's third law, we can relate rps 
to the orbital period of the secondary around the primary 
black hole Pps through: 



GMtot 2 
47r2 P= 



(18) 



where G is the gravitational constant and Mtot is the sum 
of the masses of the two black holes. 

In the observer's reference frame, the orbital period 
Ppt^ is given by: 



p — 

ps — 



(19) 



(1 + ^) 

According to iR.omero et alJ l)200rf l. density waves, 
which disturb the accretion rate and originate the superlu- 
minal components, could be produced when the secondary 
black hole crosses the primary accretion disc. As the sec- 
ondary crosses the disc twice per orbit, we used PpJ"* as 
twice the mean separation between the emergence of super- 
luminal components, which leads to Pps ~ 1.4 yr. 

Using reverberation mapping techniques, 
IPeterson et alJ il998l) estimated a virial mass for the 
central source of 3.4 x 10^ ^Q- We assumed that this value 
corresponds to the total mass of the binary system and 
from equation (18), with the values of Pps and Mtot given 
above, we found that rps ~ 5.9 x 10^^ cm. Considering 
that t he outer radius of the precessing part of t h e disc 
is rd, IPapaloizou fc TerauemI il995l) and iLarwoodl (^g^) 
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Table 5. Parameters of a possible black hole binary system in 
the inner parts of 3C 120. 

Pps (yr) rp, (cm) Mp (Mq)^ (Mq)^ Mtot (M0)<= 



1.4 



5.9 X 10^ 



3.0 X lO'^ 4.0 X 10'^ 



3.4 X 10^ 



Upper limit; 
Lower limit; 
" iPeterson et al.l <1998l) . 

the secondary, which should be higher than 4.0 x 10^ Mq. 
Summarizing, we present in Table 5 the parameters of the 
black hole binary system calculated in this section. 



Figure 4. Outer radius of the precessing disc as a function of the 
fractional mass of the primary black hole (full line). The dashed 
line indicates the separation between the primary and secondary 
black holes. 



calculated its precession period Pd in terms of the masses 
of the black holes: 



2ti 



(1- 



(t 



2n\ Gh'h 



(20) 



where n is the politropic index of the gas (e.g., n = 3/2 and 
n = 3 for the non-relativistic and relativistic cases, respec- 
tively) and 9 is the angle between the orbit of the secondary 
and the plane of the disc. 

If the jet and accretion disc are coupled, the jet pre- 
cesses at same rate than the disc {P = Pd), forming a pre- 
cession cone with half-opening angle equal to the angle of 
orbit inclination {fl = 6). Thus, replacing AIs by Mtot — Mp 
in equation (20), we can calculate rd in terms of Alp and 
Mtot: 



^d 



Stt 



\7-2nJ 



(1 + ^) 



Pcosn ^/gW, 



2/3 



1/3 



(1 



n2/3 



(21) 



where Xp = Mp/Mtot. In Fig. 4, we plot rd as a function of 
Xp using the precession parameters listed in Table 2. We can 
observe that the increase of rd with Xp is more pronounced 
after Xp ~ 0.7, such that little variations in Xp introduces 
large changes in rd. 

This formalism is valid only if the disc precesses as a 
rigid-bod y, implying that rd must be appreciably smaller 
than rps jPapaloizou fc Terauemlll995^ . This limit provides 
an additional constrain to the masses of the black holes. The 
dashed line in Fig. 4 shows the value of rps derived from 
equation (18). In order to satisfy the hypothesis of rigid- 
body precession, the allowed solutions are found for values 
of rd which lie under the dashed line in Fig. 4. It necessarily 
means that Xp < 0.88, so that Mp < 3.0 x lO'' Mq. As Mtot = 
Mp -I- Ms, we can also obtain an lower limit for the mass of 



6 CONCLUSIONS 

Based on the periodicity in the historical _B-band light curve 
and variable jet structure, we propose the existence of jet 
precession in 3C 120, with a period of 12.3 yr. 

We assume that the different apparent velocities of the 
superluminal components measured at milliarcsecond and 
larger scales are related to changes in the angle between the 
jet inlet and the line of sight due to precession, although the 
superposition of unresolved components and/or interaction 
with the environment could are acting at the largest scales. 

We show that the periodicity in the optical hght curve 
can be produced by the boosted underlying jet, with a time- 
dependent boosting factor driven by precession. An upper 
limit of 0.16 /xJy was estimated for the flux density of the 
underlying jet in the comoving reference frame. The inclu- 
sion of a secular linear term in the analysis of the long-term 
variability, as in Webb (1990), is not necessary to obtain a 
good fitting to the light curve. 

The helicoidal jet pattern found bv lWalker et alJ (1200 Jl 
is interpreted in this work also as the result of jet precession. 
The helix generated by precession reproduces quite well the 
jet aperture seen in the 1.7-GHz maps up to distances from 
the core smaller than ~80 mas, where there is a probable 
stationary component. Beyond that, the helix amplitude is 
systematically larger in the southern direction, suggesting 
the existence of an external medium that does not allow jet 
propagation. In order to produce a stationary component, 
considering a one-dimensional adiabatic relativistic jet as 
well as energy and particle flux conservation, we estimate a 
lower limit of ~12 for the ratio between jet and environment 
densities. 

Assuming that jet precession in 3C 120 is driven by a 
secondary super-massive black hole in a non-coplanar or- 
bit around the primary accretion disc, using the total mass 
of the two black holes derived from reverberation mapping 
techniques and an orbital period of approximately 1.4 yr, 
we estimate an upper limit of 3.0 x 10^ Mq for the pri- 
mary black hole mass, a lower limit of 4.0 x 10® Mq for the 
secondary mass and a separation between them of about 
5.9 X 10^^ cm. 
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APPENDIX A: OPACITY EFFECTS ON 
CORE-COMPONENT DISTANCE AND 
PRECESSION JET 



As it was pointed out prev iously (e.g.. iBlandford & KonigJ 
119791 lLobanovlir99d . Il998h . the absolute core position rcore 
depends inversely on the frequency when the core is optically 
thick, what intro duces a shift in the core-co mponent sepa- 
ration. Following i Blandford & KonigJ I^^T^, we can write 
the absolute core position as: 



4.56 X 10~^^(1 -Hz) 
D-L'y^kl^'^tp sin 

-1 2/3 



/3(1-/3C0S ?i) MP^ax/r^in) 



(Al) 



where z is the redshift, Dl is the luminosity dis- 
tance (in units of parsec), fee is a constant (fee ^ 1; 
IBlandford fc KonigJ Il97fl^ . Lsyn is the integrated syn- 
chrotron luminosity (in units of erg s~^), while Fmax and 
Fmin are related respectively to the upper and lower lim- 
its of the energy distribution of the relativistic jet particles. 
The quantities tp and v are respectively the observed aper- 
ture angle of the jet (in radians) and the frequency (in Hz); 
the former is related to the int rinsic jet aperture angle V' 
through ('e.g.. lMutel et al.lll990h : 



tan('0/2) = tan(V'72) cot 



(A2) 



The core position shift Arcoro between frequencies fi 
and 1^2 {v2 ^ i^i) is given by: 



Arcoro (i^i, 1^2) = 
Lsyn sin 



/3(l-/3cos0)ln(r,^ax/r,„in) 



4.56 X 10~^^(l + z) 
Di.'j'^kl^^ip sin (f> 

(2^2 - J^i) , 

^(maf 

U1U2 



(A3) 
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Note that if we substitute in equation (A3) a simpler 
version of eq uation (A2), ^ ~ tp' cscO, we obtain equation 
(11) given in lLobanovl l|l99Sh . 

We can see that equation (A3) depends on the angle 
between the jet and line of sight; in the case of a jet which 
is precessing, this angle is a function of time, what obvi- 
ously introduces a time dependency in Arcorc- On the other 
hand, the shifts in the core-component separations do not 
occur in a fixed direction, but they are oriented according 
to the direction whose jet inlet is pointed. As the jet inlet 
is not resolved by observations, changes in its direction will 
reflect on changes in the position angle of the core region. 
Thus, a jet component, located at a distance r from the core 
and with a position angle r), will have right ascension and 
declination offsets (Aq and A5 respectively) given by: 

AQ(//i,i/2) = r{vi)sm[r]{tahs)] - 

-Ar-core(z^i, 2^2) sin[77c(tobs)] (A4) 

/\S{ux, V2) = r{ux) cos[r;(tobs)] - 

-Arcoro(;^i, 2^2) cos[r7c(tobs)] (A5) 

where 77c is the position angle of the core in the epoch tobs in 
which observation is acquired. Known the precession model 
parameters, we are able to determine the second term of the 
equations (A4a) and (A4b) and correct the component po- 
sition by core opacity effects; using them, we can determine 
the corrected core-component distance rcorr through: 

'■corr(!^l) = ['"('^l)^ + Arcoia{vi,V2)'^ — 

-2r(i^i)Arcoro(!^i, 1^2) cos(r7 - ?7c)]^''^ (A6) 

Note that if Arcorc = 0, rcorr = r. Other particular 
case is found when there is alignment between the position 
angles of the core and of the jet component (77 = rjc), such 

that rcorr — r ^^rcore- 



APPENDIX B: GLOSSARY 

In order to facilitate the reading of the manuscript, we define 
all the symbols that appear in the text in Table Bl. 

This paper has been typeset from a TpjX/ ffl^njX file prepared 
by the author. 
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Table Bl. Definition of the symbols used in the text. 



Symbol Meaning 



G gravitational constant 

c light speed 

Hf) Hubble constant 

h //o/lOO 

qo deceleration parameter 

z redshift 

Dl luminosity distance 

r core-component distance 

Aa right ascension offset of components relative to the core position 

Ai5 declination offset of components relative to the core position 

to ejection epoch for the jet components 

ri position angle on the plane of the sky of jet components 

fi apparent proper motion of jet components 

/Sapp apparent velocity of jet components 

cf) viewing angle of jet components 

u! jet precession angular velocity 

P jet precession period 

n semi-aperture angle of the precession cone 

(po angle between the precession cone axis and the line of sight 

rjo projected angle of the precession cone axis on the plane of the sky 

P jet bulk velocity 

7 Lorcntz factor of the jet bulk motion 

<5 Dopplcr factor of jet components 

7min lower limit for the jet bulk motion 

7pB Bulk Lorentz factor of the pos-shock region 

p thermodynamic pressure 

Nj proper particle density of the jet 

Nps proper particle density of the pos-shock region 

n politropic index of the gas 

u frequency 

Sj flux density of the underlying jet in the observer's reference frame 

S' flux density of the underlying jet in the comoving reference frame 

a flux density spectral index 

Mp mass of the primary black hole 

Ms mass of the secondary black hole 

Mtot total mass inside the nuclear region 

Xp ratio between the primary and the total masses 

rps separation between the primary and secondary black holes 

Pps" orbital period of the secondary around the primary black hole in the observer's reference frame 

Pps orbital period of the secondary around the primary black hole in the source's reference frame 

outer radius of the prcccssing part of the primary accretion disc 

precession period of the accretion disk 

6 angle between the orbital plane of the secondary and the plane of the primary disc 

Lsyn integrated synchrotron luminosity 

Tmax upper limit for the energy distribution of relativistic jet particles 

Tmin lower limit for the energy distribution of relativistic jet particles 

ip' intrinsic jet aperture angle 

(/■ observed jet aperture angle 

feo constant parameter 

Arcorc angular shift in core position 

ric core position angle 



